In [1]:

    

from numpy import *
from os import path, listdir
from bokeh.plotting import figure, output_notebook, show
from myfunctions import analyze, plotS21


    

In [2]:

    

output_notebook()


    

    





    
        
        
        
    

        
        BokehJS successfully loaded.
    

In [16]:

    

datapath = '\\\\128.230.72.36\labshare\Experiments\JPM\\2QubitJPMChipDispersiveSearch\powers'
files = listdir(datapath)

filepaths = []
powers = []
results = []

for f in files:
    filepaths.append(path.join(datapath, f))
    powers.append(float(f[:3]))
    results.append(analyze(filepaths[-1]))


    

    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-54dB.s2p
f0 = 4.618029+/- 0.00001 (0.00025%) GHz
kappa = 58817424.3386+/- 167610.38068 (0.28%)
Q = 12.496+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-69dB.s2p
f0 = 4.618114+/- 0.00004 (0.00087%) GHz
kappa = 58183698.3446+/- 582998.95852 (1.0%)
Q = 12.632+/- 0.011 (0.087%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-40dB.s2p
f0 = 4.617991+/- 0.00001 (0.0002%) GHz
kappa = 58713262.7773+/- 130840.12286 (0.22%)
Q = 12.518+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-59dB.s2p
f0 = 4.617989+/- 0.00001 (0.00029%) GHz
kappa = 59474564.29+/- 192825.84690 (0.32%)
Q = 12.358+/- 0.011 (0.089%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-70dB.s2p
f0 = 4.618105+/- 0.00004 (0.00089%) GHz
kappa = 59286504.2344+/- 598837.42662 (1.0%)
Q = 12.397+/- 0.011 (0.089%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-58dB.s2p
f0 = 4.617986+/- 0.00001 (0.00029%) GHz
kappa = 59389419.8117+/- 194522.60492 (0.33%)
Q = 12.376+/- 0.011 (0.089%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-52dB.s2p
f0 = 4.618008+/- 0.00001 (0.00024%) GHz
kappa = 59035843.6069+/- 158351.26921 (0.27%)
Q = 12.450+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-65dB.s2p
f0 = 4.618044+/- 0.00003 (0.00064%) GHz
kappa = 59563011.4111+/- 428890.58705 (0.72%)
Q = 12.340+/- 0.011 (0.089%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-56dB.s2p
f0 = 4.618038+/- 0.00001 (0.00025%) GHz
kappa = 58808391.3839+/- 168702.72811 (0.29%)
Q = 12.498+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-61dB.s2p
f0 = 4.61804+/- 0.00002 (0.00034%) GHz
kappa = 59080253.0406+/- 227479.10628 (0.39%)
Q = 12.440+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-20dB.s2p
f0 = 4.617975+/- 0.00001 (0.00019%) GHz
kappa = 58426850.5841+/- 126562.27355 (0.22%)
Q = 12.579+/- 0.011 (0.087%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-67dB.s2p
f0 = 4.618089+/- 0.00004 (0.00081%) GHz
kappa = 59294895.4232+/- 539697.71970 (0.91%)
Q = 12.396+/- 0.011 (0.089%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-53dB.s2p
f0 = 4.618008+/- 0.00001 (0.00023%) GHz
kappa = 58881107.4434+/- 151110.30419 (0.26%)
Q = 12.482+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-68dB.s2p
f0 = 4.618059+/- 0.00004 (0.00086%) GHz
kappa = 59399094.0853+/- 576913.91448 (0.97%)
Q = 12.374+/- 0.011 (0.089%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-57dB.s2p
f0 = 4.618021+/- 0.00001 (0.00032%) GHz
kappa = 59120456.1218+/- 211568.04167 (0.36%)
Q = 12.432+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-51dB.s2p
f0 = 4.618015+/- 0.00001 (0.00022%) GHz
kappa = 58881640.9702+/- 146202.69503 (0.25%)
Q = 12.482+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-66dB.s2p
f0 = 4.61797+/- 0.00003 (0.00071%) GHz
kappa = 58994131.3515+/- 475866.24280 (0.81%)
Q = 12.458+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-64dB.s2p
f0 = 4.61804+/- 0.00003 (0.00057%) GHz
kappa = 59477299.5844+/- 383992.31857 (0.65%)
Q = 12.357+/- 0.011 (0.089%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-55dB.s2p
f0 = 4.618022+/- 0.00001 (0.00025%) GHz
kappa = 58621558.8079+/- 164353.58327 (0.28%)
Q = 12.538+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-63dB.s2p
f0 = 4.618032+/- 0.00003 (0.00056%) GHz
kappa = 58922556.9749+/- 375348.39174 (0.64%)
Q = 12.474+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-60dB.s2p
f0 = 4.61801+/- 0.00001 (0.0003%) GHz
kappa = 58970966.0293+/- 198248.21678 (0.34%)
Q = 12.463+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-62dB.s2p
f0 = 4.618034+/- 0.00002 (0.00045%) GHz
kappa = 58788996.8444+/- 300254.95294 (0.51%)
Q = 12.502+/- 0.011 (0.088%)






    















    




\\128.230.72.36\labshare\Experiments\JPM\2QubitJPMChipDispersiveSearch\powers\-50dB.s2p
f0 = 4.617991+/- 0.00001 (0.00022%) GHz
kappa = 58994436.9126+/- 147380.55750 (0.25%)
Q = 12.458+/- 0.011 (0.088%)






    

In [35]:

    

x = []
y = []
for i, result in enumerate(results):
#     print(result)
    x.append(powers[i])
    y.append(results[i]['f0']*1e-9)


    

In [25]:

    

from bokeh.plotting import Plot


    

In [50]:

    

p = figure(title="Resonance at various Powers", x_axis_label="Resonance (GHz)", y_axis_label="Power (dB)")
p.circle(y,x)
show(p)


    

In [52]:

    

p = figure(title="Tiny dispersive shift, high and low powers", x_axis_label="Frequency (GHz)",y_axis_label="S21 (dB)");
for i, result in enumerate(results):
    if powers[i] == -20 or powers[i] == -70:
#         print(powers[i])
#         print(results[i])
        p.line(result['freq']/1e9, result['dBm'], legend = str(powers[i]) + ' dB')
show(p)


    

In [ ]: